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ABSTRACT 

Querying the Spitzer archive for optically-selected active galactic nuclei (AGN) observed in high-resolution-mode spec- 
troscopy, we identified radio and/or interacting galaxies with highly turbulent motions of Eb gas at a temperature of 
a few hundred Kelvin. Unlike all other AGN that have unresolved H2 line profiles at a spectral resolution of ~600, 
3C236, 3C293, IRAS09039+0503, MCG- 2-58-22 and Mrk463E have intrinsic velocity dispersions exceeding 200 kms" 1 
for at least two of the rotational SO, SI, S2, and S3 lines. In a sixth source, 4C12.50, a blue wing was detected in the 
SI and S2 line profiles, indicating the presence of a warm molecular gas component moving at -640 kms -1 with respect 
to the bulk of the gas at systemic velocity. Its mass is 5.2 xlO 7 Mq, accounting for more than one fourth of the H2 
gas at 374K, but less than 1% of the cold H2 gas computed from CO observations. Because no diffuse gas component 
of 4C12.50 has been observed to date to be moving at more than 250 kms -1 from systemic velocity, the H2 line wings 
are unlikely to be tracing gas in shock regions along the tidal tails of this merging system. They can instead be tracing 
gas driven by a jet or entrained by a nuclear outflow, which is known to emerge from the west nucleus of 4C12.50. 
It is improbable that such an outflow, with an estimated mass loss rate of 130 Mqjx ~ , entirely quenches the star 
formation around this nucleus. 
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1. Introduction 

Large-scale feedback effects such as jets and outflows from 
active galactic nuclei (AGN) are thought to be capable of af- 
fecting the formation of new stars in their host galaxies. The 
triggering of star formation by compression of gas (e.g., van 
Breugel et al., 1985), as well as the suppression of star for- 
mation by heating of gas that prevents its further collapse 
(e.g., Nesvadba et al., 2010) have been observed in local 
AGN. Cosmological simulations have suggested that AGN 
feedback effects, which are often associated with mergers, 
could make galaxies appear red, or even explain the ob- 
served luminosity functions of galaxies (Croton et al., 2006; 
Hopkins et al. , 2006) . Combined with multi- wavelength ob- 
servations indicating that the star-formation history and 
the black-hole-accretion history of the Universe peak at 
comparable redshifts, between l<z<3 (Marconi et al., 2004; 
Merloni et al., 2004), this suggests that AGN feedback could 
have affected the shape of present-day galaxies consider- 
ably. 

Extensive tests of the role of AGN feedback on the in- 
terstellar medium (ISM) of galaxies require a detailed kine- 
matic study of outflowing gas in local sources. Signs of mas- 
sive gas outflows have been detected for ionized atomic gas 
(e.g., Veilleux et al., 1995; Emonts et al., 2005; Holt et al., 
2006; Muller-Sanchez et al., 2006), neutral atomic gas (e.g., 
Morganti et al., 2005; Rupke et al., 2005), and molecular 
CO and OH gas (e.g., Curran et al., 1999; Das et al., 2005; 
Garcfa-Burillo et al., 2009; Sakamoto et al., 2009; Feruglio 



et al., 2010; Fischer et al., 2010; Sturm et al., 2011). In this 
letter we present evidence for the first detection of highly 
turbulent motions of H2 gas at a temperature of a few 
hundred Kelvin as seen with Spitzer for several local AGN. 
We adopt H o =70 kms" 1 Mpc" 1 , O M =0.3, and Q A =0.7 
throughout this work. 

2. The sample selection 

We queried for turbulence in the warm H2 gas in local 
AGN using mid-IR spectra obtained with Spitzer in high- 
resolution mode. The full archival sample comprises 298 
sources that are classified as AGN based on optical spec- 
troscopy. It is presented in Dasyra et al. (2011) together 
with the data reduction techniques and the extracted spec- 
tra. 

To look for turbulent H2 gas motions, we examined the 
profiles of the purely rotational (0-0)S0 28.22 rim, (0-0)Sl 
17.04 /j,m, (0-0)S2 12.28 fun, and (0-0)S3 9.66 11m lines. We 
searched for either resolved profiles with velocity disper- 
sion (7>200 kms -1 , or for profiles with asymmetric wings 
that are characteristic of outflows. To ensure the reliability 
of our results, we only examined sources with at least two 
lines of signal-to-noise (S/N) ratio >5. We also requested 
that at least two lines suggest a similar kinematic pattern, 
i.e. a wing or a resolved profile. To consider a line resolved 
we requested that its full width at half maximum (FWHM) 
value minus the FWHM error exceeds the instrumental res- 
olution R plus the resolution error at the observed-frame 
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Fig. 1. Normalized, continuum-subtracted molecular and ionized gas line profiles in local AGN with resolved H2 emission. The 
filled gray area represents a Gaussian with FHWM equal to the resolution of the IRS. 



Table 1. Fluxes and widths of resolved H2 lines in local AGN. 



Source 


z 


fsi 


FWHMsi" 


fs2 


FWHMsj" 


f S3 


FWHMs3° 




T 


M cold " 






HT 17 W m~ 2 


kms -1 


HT 17 W m- 2 


kms -1 


HT 17 W m~ 2 


kms -1 


10 7 M s 


K 


10 9 M 


3C236 


0.0989 


1.06±0.07 


582 (750) 


0.49±0.08 


557 (719) 


0.81±0.09 


624 (767) 


6.10 


345 


<5.1 


3C293 


0.0450 


5.30±0.38 


510 (716) 


1.79±0.31 


553 (732) 


2.96±0.38 


519 (714) 


6.62 


323 


23 


IRAS09039+0503 


0.1254 


3.09±0.54 


519 (700) 


1.34±0.06 


575 (730) 


2.04±0.13 


569 (751) 


31.6 


334 




MCG-2-58-22 


0.0472 


3.86±0.42 


544 (740) 


1.50±0.25 


636 (796) 
... (...) 


3.40±0.46 


662 (824) 


4.19 


359 


5.7 


Mrk463E 


0.0507 


3.37±0.59 


556 (748) 


<1.98 


3.76±0.44 


590 (767) 


4.06 


378 


1.2 


4C12.50 (main c ) 


0.1217 


1.64±0.21 


... (519) 


1.20±0.11 


... (552) 


1.78±0.46 


... (652) 


13.9 


374 


15 


4C12.50 (wing c ) 


0.1196 ti 


0.62±0.21 


... (568) 


0.44±0.11 


521 (690) 




... (...) 


5.19 


374<= 





Notes. None of these sources had a reliable SO detection owing to the high noise levels in the long wavelength array. Rest-frame, 
instrumental-broadening-corrected FWHM. Their measured values appear in parenthesis. Their error bars correspond to 5—10% 
of their values. 6 Cold H2 gas masses calculated from CO observations (Evans et al., 1999, 2002; Bertram et al., 2007; Saripalli & 
Mack, 2007) using an H 2 / CO mass conversion factor of 1.5M /(K kms _1 pc 2 ) for the IR-bright systems Mrk463E and 4C12.50 
(Evans et al., 2002), and a standard Galactic value of 4.8M0/(K kms~ 1 pc 2 ) for all other sources. They are converted to the 
adopted cosmology. The main and wing components correspond to the primary and the secondary Gaussian functions in 

Figure 2. ^ The Gaussians that best fit the SI and S2 line wings peak at -646 kms -1 and -634 kms -1 from systemic velocity, 
respectively. ' e ' The excitation temperature of the rapidly moving H2 gas is assumed to be identical to that of the H 2 at systemic 
velocity to facilitate the comparison of their masses. 



wavelength of the line (Dasyra et al., 2011). The average R 
value in the 12.0— 18.0 ^tm range, which comprises the H2 
SI and S2 transitions, is 507±66 kms -1 . 

3. Results: Sources with highly turbulent or rapidly 
moving II L > gas in the warm phase 

Of the 298 sources 62 had at least two H2 lines de- 
tected with S/N>5. The profiles of two or more lines were 
spectrally resolved in only five sources, namely 3C236, 
3C293, IRAS09039+0503, MCG-2-58-22, and Mrk463E 
(Fig. 1). Their velocity dispersions are in the range 
220<(7<280kms~ 1 (Table 1). Because all these sources are 
radio galaxies and/or interacting systems, the high turbu- 
lence in their warm H2 motions can be driven by AGN 
feedback mechanisms, by gravitational instabilities, or by 
supernova winds. Still, no mechanism is efficient enough to 
kincmatically distort a detectable mass of warm H2 gas 
to velocity dispersions exceeding 300 kms -1 . We computed 
the excitation temperature T and mass of the turbulent gas 
(Table 1) using the detected SI, S2, and S3 line fluxes as 
in Higdon et al. (2006). At temperatures of 300-400 K, its 
mass is typically on the order of 1% of the cold H2 gas mass 
indicated by CO observations. 

Further outflow or inflow signatures were sought for 
in the H 2 line wings and in the difference of the H 2 
recession velocity from the systemic velocity, V sys . The 
latter was determined from the [Nell] line, emitted by 



ions that are abundant in star-forming regions and in 
the AGN vicinity owing to their low ionization poten- 
tial, 21.56 eV. The H2 recession velocity agreed within 
the errors with V sys for all sources, including those 
with massive outflows of the gas that is photoionized by 
the AGN and that is traced by the [NeV] or [OIV] 
lines (i.e., 3C273, IRAS13342+3932, IRAS05189-2524, 
IRAS15001+1433, IRAS23060+0505, Mrk609; Dasyra et 
al., 2011). 

The only source with wing signatures in its H2 line 
profiles was 4C12.50, also known as IRAS13451+1232 or 
PKS1341+12. The SI and S2 line wings (Fig. 2; upper 
panel), detected with S/N>3, point at two warm molec- 
ular gas kinematic components in this source. The peak of 
the secondary Gaussian that is needed to fit both profiles is 
blueshiftcd by ^640 kms -1 from the primary Gaussian at 
V sys . The flux ascribed to the primary Gaussian is only 
2.6 and 2.7 times higher than that ascribed to the sec- 
ondary Gaussian of the SI and S2 transitions, respectively. 
Using the actual flux in each Gaussian (Table 1), we sepa- 
rately calculated the H2 gas mass for each kinematic com- 
ponent. We find the excitation temperature of the bulk of 
the gas to be 374±12 K. It is equal to the inverse of the 
slope of the line that best fits the excitation diagram points 
(Rigopoulou et al., 2002). The excitation diagram (Fig. 2; 
lower panel) shows the natural logarithm of the number of 
electrons n e that descended from the upper to the lower 
state, normalized by the statistical weight g of the transi- 
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Fig. 2. Upper panel: Normalized, continuum-subtracted molec- 
ular and ionized gas line profiles for 4C12.50. In addition to 
the profiles of the MIR lines, the 5007 A [O III] profile that is 
convolved to the resolution of IRS is presented for comparison. 
The blue wing is detected in all lines except for the unresolved 
S3. Lower panel: H2 excitation diagram constructed separately 
for the primary Gaussian component (open diamonds) and the 
secondary Gaussian component (stars) of each H2 line. 



tion, as a function of the temperature that corresponds to 
the energy of the upper state E up divided by the Boltzmann 
constant k. The value of n e is computed as L/ (ahv), where 
h is the Planck constant, a is the Einstein coefficient of 
the transition, and v is the frequency of the emitted line. 
For a single temperature of 374K we find that the mass of 
the H2 gas at systemic velocity is 1.4 x10 s M (see also 
Higdon et ah, 2006). Assuming (for simplicity) that T is the 
same for both H2 kinematic components (see Figure 2), we 
find that the mass of the H2 gas moving at 640 km s -1 is 
5.2xl0 7 M . This is 0.3% of the cold H 2 gas mass in the 
west nucleus of 4C12.50, which was found from CO obser- 
vations to be 1.5xlO lo M (Evans et al., 2002). 

4. Discussion: An AGN-driven molecular gas 
outflow in 4C12.50? 

Even though 4C12.50 is an IR-bright system of two inter- 
acting galaxies (Axon et al., 2000), the secondary Gaussian 
of Figure 2 is unlikely to be tracing gas in the east nucleus, 
which is mostly located outside the IRS slit (Fig. 3). Any 
residual gas from the east nucleus inside the slit would be 
moving at a velocity comparable to the difference in the 
recession velocity of the two nuclei, ~200kms~ 1 . This dif- 
ference is computed from ionized gas kinematics (Holt et 
ah, 2003), and it is confirmed by stellar kinematics from 
CO absorption features presented in Dasyra et al. (2006). 
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Fig. 3. Composite Ha (cyan) and 5900 A continuum (red) im- 
age of 4C12.50, constructed from Hubble Space Telescope data 
(Batcheldor et al., 2007). To enhance the visibility of low surface 
brightness structures, we removed all symmetric galaxy compo- 
nents (i.e., a common bulge and two residual disks located at the 
position of the two nuclei) from the 5900 A image with GALFIT 
(Peng et al., 2002). Both disks were found to have an inclina- 
tion of ~20°. Irregularly moving components are found within 
the IRS slit, marked with a yellow box for both nod positions. 
(A) marks the position of two super star clusters with veloci- 
ties that are blueshifted by up to ^250 kms -1 from systemic 
velocity (Rodriguez Zaurin et al., 2007). 



Besides the two nuclei residing in a common bulge, 4C12.50 
also has off-nuclear gas concentrations and super star clus- 
ters in tidal tails. The blueshifted H2 emission could arise 
from gas in a tidal tail inside the IRS slit (Fig. 3), which 
is entering a shock front created during the collision (see 
Cluver et ah, 2010). The tidal tail could be moving faster 
or be at a different inclination angle i than its corresponding 
nucleus, of ^20° in either case. Because the deprojected tail 
velocity would be equal to 640 / sin(i) kms -1 , it could reach 
a value as high as ^2000 kms -1 . This scenario is also un- 
likely given that no off-nuclear, large-scale (<20 kpc) kine- 
matic component has been observed to be moving faster 
than ±250 km s -1 from V sys along the line of sight (Holt 
et ah, 2003; Rodriguez Zaurin et ah, 2007). 

A scenario that agrees better with existing observations 
is that the H2 gas is moving toward us driven by feed- 
back mechanisms (e.g., Alatalo et ah, 2011). Optical spec- 
troscopy indicated the presence of an outflow from the west 
nucleus of 4C12.50 by revealing the existence of three kine- 
matic components for the nuclear [OI], [S II], and [OIII] 
emission (Holt et ah, 2003). Most of the [OIII] emission is 
blueshifted by 400 kms -1 , while its broadest component (of 
FWHM-1900 kms- 1 ) is blueshifted by 2000 kms" 1 from 
V sys . MIR spectroscopy suggested an AGN- related nuclear 
outflow of ionized gas. Blue wings were observed in the 
profiles of the [OIV] 25.89 fim and the [NeV] 14.32 //m 
lines (Spoon & Holt, 2009; Dasyra et ah, 2011), emitted by 
ions that are primarily found in hard radiation fields. Radio 
observations revealed an HI absorption line blueshifted by 
~1000 kms -1 from V sys (Morganti et ah, 2004). Because 
a background radio source is required for HI absorption 
to be seen, the hydrogen clouds are likely to be located 
between the observer and the AGN or its jet. Estimates 
of the outflow mass range from 8xlO 5 M for the ionized 
10 4 K gas (Holt et ah, 2011) to 5.6xlO 8 M for the neu- 
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tral gas traced by NalD (Rupke et al., 2005), bracketing 
our mass estimate, M out , of 5.2 xlO 7 M Q for the outflowing 
-400K H 2 gas. 

If the outflow is caused by AGN radiation pressure or 
winds (Holt et al., 2011), it can be considered spherical. 
If the gas is also distributed in a sphere, and if its den- 
sity is falling with the inverse square of the distance from 
the center, its mass outflow rate, M, will be given from 
the product M out V^ U iii _1 . For an outflow velocity V out of 
640 km s -1 and for a radius R of 270 pc, as estimated for 
the CO gas assuming that it is thermalized to the dust 
temperature (Evans et al., 2002) and as converted to the 
adopted cosmological distance, M will be 130 M Q yr _1 . An 
outflow of these properties is unlikely to entirely suppress 
star formation in the west nucleus of 4C12.50, whose star- 
formation rate (SFR) is estimated to be between 370-1380 
M Q yr . The lower value is found from the CO mass using 
a gas consumption timescale of 4xl0 7 yrs, while the up- 
per value is found by folding the CO-based H 2 mass and 
radial extent in the Schmidt-Kennicutt law (Evans et al., 
2002). If the outflow were symmetric, both a blue and a 
red H2 line profile wing should exist unless the gas moving 
away from the observer is obscured by dust. This is plau- 
sible for a source with E(B-V) of 1.44 magnitudes (Holt et 
al., 2011) and 9.7/im optical depth of 0.59 (Veilleux et al., 
2009), which could translate into a 10—20% absorption of 
the total flux at 17 /im (Li & Draine, 2001), and which could 
preferentially suppress the red line wing for a circumnuclear 
dust distribution. 

Alternatively, a radio jet encountering clouds on its 
path could be driving the outflow (e.g., Dietrich & Wagner, 
1998). A jet is indeed known to emerge from the west nu- 
cleus of 4C12.50. It extends out to 45 pc in projection in the 
north and 170 pc in the south (Stanghellini et al., 1997), 
and it propagates close to the speed of light at a small an- 
gle from the line of sight (Lister et al., 2003). A previous 
flare of the jet, undetected in the radio because of its weak 
signal, could be associated with the shocked gas whose ex- 
tended X-ray emission peaks at 20 kpc south of the nucleus 
(Sicmiginowska et al., 2008). The scenario of a jet-driven 
outflow can easily explain the observed line profiles. The 
detection of a blue or a red wing is random since it de- 
pends on the location of the clouds with respect to the jet 
propagation axis. 



5. Summary and concluding remarks 

We queried the archival catalog of 298 optically-selected 
AGN observed with Spitzer IRS in high-resolution mode 
(Dasyra et al., 2011), aiming to identify sources with tur- 
bulent motions of their warm molecular gas. We examined 
the profiles of the H 2 SO, SI, S2, and S3 lines and found 
only five radio and/or interacting galaxies with <r>200 
kms -1 but no source with <r>300 kins" 1 . In a sixth source, 
4C12.50, the SI and S2 lines have a blue wing that points 
at warm gas moving toward us with 640 kms -1 . Its mass, 
5.2 xlO 7 Mq, corresponds to an impressively high fraction, 
^1/4, of the total ~400K H 2 gas mass. While it could be 
tracing shock regions along tidal tails, it is more likely to 
be tracing an AGN jet or wind-driven outflow, known from 
ionized and neutral gas kinematic studies. Even if all of 
this gas is entrained by an outflow, it is unlikely to entirely 
suppress star formation in 4C12.50. Additional tests of the 



role of AGN feedback mechanisms in increasing the tur- 
bulence of the molecular gas require observations of high- 
rotational-number transitions of CO molecules that can be 
mostly excited by the AGN (van der Werf et al., 2010). An 
essential role in revisiting this question will also be played 
by the Atacama Large Millimeter Array, which, via high- 
resolution studies of the cold gas, will enable a comparison 
never performed before: the computation of the warm-to- 
cold molecular gas mass ratio in an outflow vs the rest of 
the ISM. 
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